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Abstract: We propose a duality between quiver gauge theories and the combinatorics of 
dimer models. The connection is via toric diagrams together with multiplicities associated 
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construction of the toric space). These multiplicities may be computed from both sides and 
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arbitrary orbifolds of a toric space, and allow a new algorithmic method for exploring the 
phase structure of the quiver gauge theory. 
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1. Introduction 

Dimer models 1 are two-dimensional statistical mechanical models of tilings of graphs. They 
have been recently introduced to the string theory literature in [IJ, which proposed to relate 
the dimer models to the target space topological string held theory. This viewpoint has been 
developed |2[ for the crystal models also discussed in [|I|, but the precise topological string 
theory interpretation of the dimers has not yet been made clear 2 . 

In this paper we develop a different string theoretical interpretation of the dimer models, 
and relate them to the gauge theories that live on the world-volume of D-branes probing non¬ 
compact toric spaces. These theories are generally known under the name of quiver gauge 
theories and were the subject of study of many recent papers. A generic feature of these 

1 dimer (n): a compound whose molecules are composed of two identical monomers. 

2 See however the recent papers || Q] which pointed out that the Ronkin function, which is associated 

to the limit curve of the dimer model, can be used to recover the instanton contributions to the topological 
string partition function. 
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theories is that for a given singular non-compact manifold there exist infinitely many quiver 
gauge theories that probe the same singularity. These theories are all related by Seiberg 
duality and form a duality tree that encodes the relations between theories. Examples of 
such trees and a discussion of this phenomenon can be found in [|J. Out of all the theories 
in the duality tree there is a special subclass of theories which are called the toric phases 
of the quiver gauge theory. In all known examples this subclass contains a finite number of 
theories. The toric phases are characterized by having the same rank for all gauge groups in 
the quiver, with each field appearing exactly twice in the superpotential. A discussion of all 
of the toric phases for a class of quiver gauge theories called Y p,q was recently performed in 

i- 

This paper is based on the observation that the coefficients appearing in the determinant 
of the Kasteleyn matrix of the graph associated to a dimer model are the same numbers that 
count the multiplicities of fields in the linear sigma model realization of the toric variety. 
Moreover, the various toric phases of the quiver gauge theories, which are related by toric 
duality, may all be obtained from related dimer models. Thus, the dimer models know about 
the phase structure of the gauge theories living on the moduli space of D-branes probing 
toric singularities. This paper explores some aspects of this correspondence. 

In section || we revise some basics of quiver gauge theories, and the forward and inverse 
algorithms for going between the quiver theory and its’ moduli space. In section [3] the dimer 
models are introduced and their relation to toric geometry is elucidated. In terms of the 
dimer models, certain aspects of the toric geometry become clear: for example they provide 
simple formulae for the linear sigma model field multiplicities of orbifolds of any non-compact 
toric 3-fold, which recover the known examples in particular cases. 

We show how the operation of Higgsing the gauge theory to remove points from the toric 
diagram 3 corresponds to cutting edges of the graph, and describe how this procedure may 
be used to systematically enumerate the toric phases of a given singularity. 

In a subsequent paper we will explore the map between the quiver gauge theory and 
dimer pictures and combine it with the results from the dimer side to obtain a deeper 
understanding of the quiver gauge theories and their embedding in string theory. 

2. Review of quivers and toric singularities 

When D-branes sit on a non-compact singular Calabi-Yau (CY) manifold, X, the world vol¬ 
ume theory living on them is a quiver gauge theory which involves gauge fields, matter fields 
and interaction terms in a theory that is generically supersymmetric with 4 supercharges. In 
special cases when the manifold has more symmetries, the supersymmetry may be enhanced 
to as many as 16 supercharges. Taking the number of D-branes to be N, examples are 

3 also called “partial resolution” in previous string theory literature. 
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X = C 3 with quiver gauge theory SYM U(N) and 16 supercharges; X = C x C 2 /Z n with 
quiver gauge theory having 8 supercharges and specihed by the affine Dynkin diagram of 
A n -h etc. 

The problem of finding a quiver gauge theory for a given toric singularity is a difficult 
problem and is still not fully understood. One of the methods in constructing this quiver 
gauge theory is known under the name “partial resolution of Abelian orbifolds” Bi- An 
algorithm for computing the gauge theory for a given toric singularity was used in |9j to find 
the gauge theories corresponding to all of the toric del Pezzo surfaces. 4 

The sketch of the algorithm goes as follows. To every toric diagram one can embed it 
in a sufficiently large toric diagram with a known quiver gauge theory. Such a diagram can 
be chosen to be C 3 /Z n x 7L m with n and m the smallest integers such that the orbifold toric 
diagram (which looks like a triangle with sides of dimension (n+ 1) x (m + 1)), contains the 
toric diagram of interest. Then one can remove points by introducing FI terms to the cor¬ 
responding gauge theories. This leads to a Higgs mechanism and generation of mass terms 
for some of the fields, which can then be integrated out. As we follow through the process 
of removing points one can end up with the right quiver and superpotential for the desired 
toric diagram. While this is a simple process to describe, it becomes impractical to imple¬ 
ment for large toric diagrams. With dimers, however, the amount of computation reduces 
substantially and it is possible to proceed with this algorithm for larger toric diagrams. See 
section |5] for more details. 

One interesting aspect of the algorithm is its non-uniqueness in producing a quiver gauge 
theory, thus leading to a one-to-many map between a geometric singularity and quiver gauge 
theories. The corresponding quiver gauge theories are different UV theories but flow to the 
same universality class in the IR, characterized by the singular geometry. The non-uniqueness 
phenomenon is called “Toric duality” U and is found to be equivalent to Seiberg duality 
m 0 0- Each different quiver theory resulting from the algorithm is called a “toric 
phase” and has the property that all ranks of the gauge groups are equal to the number of 
D-branes, N, and every bi-fundamental held appears in the superpotential precisely twice: 
once with a positive sign and once with a negative sign. 

The computation of the gauge groups and matter fields turns out to be a simpler task 
than that of computing the superpotential. Indeed a few methods, different and even faster 
than the inverse algorithm of [|9j, were used in the literature to compute the matter content 
of the quiver theory: using (p,q) webs JO], Ext groups |T^|, exceptional collections 
and others. Recently a generic computation of superpotentials was proposed in [[f~5|] . One 
of the results of the present work in relating dimers to quiver gauge theories is a major 

4 In the literature the process of finding the toric data for a given quiver gauge theory is called the 
“Forward Algorithm” and the process of finding quiver gauge theories from a given toric data is called the 
“Inverse Algorithm”. 
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simplification in computing superpotentials. If, using independent data, one comes up with 
a guess for a quiver gauge theory, then one of the immediate subsequent computations to 
be performed is the construction of its’ moduli space of vacua. This moduli space coincides 
with the CY manifold probed by the D-branes. The dimer methods that we introduce in 
this paper allow a quick computation of the toric data of this CY manifold, much faster than 
the forward algorithm or any other existing methods do. 

Another interesting aspect of the algorithm used in || is the appearance of a multi¬ 
plicity of fields in the gauged linear sigma model. A detailed study of these multiplicities 
and of toric duality, done in |16|, proposes that the different toric phases of a given singular 
manifold are characterized by different multiplicities. All known examples studied so far 
have the property that the multiplicities of the toric diagram are different for different toric 
phases. Using the methods presented in this paper it becomes clear that there are coun¬ 
terexamples to this proposal as the results on Y 6,0 indicate. Furthermore these multiplicities 
have beautiful combinatorial properties and the explicit formulas given in this paper shed 
light on their origin. 

The “Forward Algorithm” starts with a quiver gauge theory for one D-brane, (N = 1). 
In the language of 3+1 dimensional supersymmetry with 4 supercharges it has n U{ 1) vector 
multiplets, m bi-fundamental chiral multiplets and m — n terms in the superpotential. At 
the end of the computation one ends up with c fields Pi 5 subject to c — 3 moment map 
relations. This gives a (c - 3) x c charge matrix which is called Qt■ The integral co-kernel 
of Q t is a 3 x c matrix called G t with columns corresponding (up to repetition) to nodes of 
the three dimensional toric diagram. All these nodes lie on a 2 dimensional sub-lattice, due 
to the CY condition. Due to the repetition, each node comes with a multiplicity and it is 
these multiplicities that form the main interest of the present paper. Further details on the 
forward algorithm can be found in p, [TEj and in [|TU] . 

The forward algorithm uses linear algebra and integer programing and turns out to be 
a useful tool in computing toric data, including their multiplicities, for sufficiently small 
quivers. Here small means small n and small m. Once the size of m and n become large it 
becomes exponentially difficult to compute the multiplicities and hence practically impossi¬ 
ble. For example the computation of the multiplicities for the quiver theory of C 3 /Z 4 x Z 4 
which has n = 16 and m = 48 would take the order of few days on an average-speed desktop 
computer. Obviously it becomes clear that we need more efficient methods to compute these 
multiplicities for generic toric diagrams. As we will see in this paper, dimers provide the right 
tool to compute these multiplicities with arbitrary length and with analytic expressions. In 
particular for orbifolds of C 3 , of the conifold, or of any other type of singularity, we have 
explicit expressions which were found by Kenyon [|20| and are given in equation ( [3.131 ) and 


5 In the math literature the fields pt are the homogeneous coordinates of the toric variety. See |17 
detailed discussion on the C 3 /Z 2 x Z 2 orbifold. 
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generalized in ( 3.14 ). 

The collection of multiplicities which were known so far were summarized in |1(| . Some 
others were found in |^T|. A few more, found using the forward algorithm, will be published 
in an upcoming work [19]. One of the points of this paper is the ability to generalize all 


these results into a simple computational tool which produces many more new results. 

In a recent paper |22| the quiver gauge theories of the Y p,q manifolds were found. A 
subsequent work [J6J computed all the toric phases for each Y p,q model and demonstrated the 
rich structure they have. One non-trivial check for the duality between the Y p,q manifolds 
and their quiver gauge theories was performed in computing volumes of supersymmetric 
cycles and comparing them to i?-charges computed in the gauge theory. Another check 
which follows from this work is to perform the computation of the moduli space of vacua of 
the gauge theories using their dimer description. A subsequent work |23| will present another 
infinite family of quiver gauge theories which have no known metrics. The methods in the 
current paper will provide a consistency check for the toric description of the corresponding 
CY manifolds. 

The geometrical meaning of these multiplicities still remains unclear, and it remains 
unknown whether they are connected in any way to topological strings. We expect that 
studying such connections will be a highly fruitful area of research. 

This work introduces a connection between dimers and multiplicities of toric diagrams, 
relevant to quiver gauge theories. In a subsequent publication 0 we will explore the 
relationship between dimers and quiver gauge theories in greater detail, by extending the 
correspondence between quiver gauge theories and toric geometries as their moduli space. 
Dimers turn out to be the right tool to study both in a simple and efficient fashion. 


3. Dimer models 

This section is based on [|(], [25j. 

A bipartite graph is a graph with the property that all vertices can be coloured black 
or white, such that every black vertex is only adjacent to white vertices, and vice versa. A 
perfect matching of a bipartite graph is a subset of edges (“dimers”) such that every vertex 
in the graph is an endpoint of precisely one such edge. A dimer model is the statistical 
mechanics of such a system, i.e. of random perfect matchings of the graph with assigned 
edge weights. We will be interested in dimers on doubly-periodic graphs, i.e. defined on the 
torus T' 2 . 

Many important properties of the dimer model are governed by the Kasteleyn matrix 
K(z , w), a weighted, signed adjacency matrix of the graph with (in our conventions) the rows 
indexed by the white nodes, and the columns indexed by the black nodes. It is constructed 
as follows: 


5 
















To each edge in the graph, multiply the edge weight by ±1 so that around every face 
of the graph the product of the edge weights over edges bounding the face has the following 
sign 


sign (I] 


ed = 


+ 1 if (# edges) = 2 mod 4 
— 1 if (# edges) = 0 mod 4 


(3.1) 


It is always possible to arrange this [^ 6 [ . 

The colouring of vertices in the graph induces an orientation to the edges, for example 
the orientation “black” to “white”. Construct paths 7 ™, in the dual graph that wind 
once around the (0,1) and (1, 0) cycles of the torus, respectively. For every edge crossed by 
7 , multiply the edge weight by a factor of w or 1/w (respectively z, 1 /z) according to the 
relative orientation of the edges in G crossed by 7 . 

The adjacency matrix of the graph T weighted by the above factors is the Kasteleyn 
matrix K(z,w) of the graph. The determinant of this matrix P(z,w ) = det K is called the 
characteristic polynomial of the graph, and will be our main object of study. 

Given an arbitrary reference matching M 0 on the graph, the difference M — M 0 defines 
a set of closed curves on the graph in T 2 . This in turn defines a height function on the faces 
of the graph: when a path in the dual graph crosses the curve, the height is increased or 
decreased by 1 according to the orientation of the crossing. A different choice of reference 
matching M 0 shifts the height function by a constant. Thus, only differences in height are 
physically significant. 

In terms of the height function, the characteristic polynomial takes the following form: 


P(z,w) = z' 


- yjhyO 




c h x ,hy{ 1 ) 


hx+hy hx hy 


(3.2) 


where Ch xi h y are integer coefficients that count the number of paths on the graph with height 
change (h x ,h y ) around the two fundamental cycles of the torus. These coefficients do not 
have an obvious geometrical meaning in toric geometry 6 , but will turn out to be equal to the 
multiplicity of holds required to engineer the toric variety as the moduli space of a D-brane 
world-volume quiver gauge theory. 

In the following it will be convenient to define the following quantity: 


c = 



(3.3) 


which counts the total multiplicity of the linear sigma model fields. This number was com¬ 
puted for special cases 0 , 0 . 0 . 

The overall normalization of P(z, w) is not physically meaningful: since the graph does 
not come with a prescribed embedding into the torus (only a choice of periodicity), the paths 

6 There should be an interpretation in terms of the symplectic quotient construction of the space. 
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lz,w winding around the primitive cycles of the torus may be taken to cross any edges en 
route. Different choices of paths 7 multiply the characteristic polynomial by an overall power 
z l wi, and by an appropriate choice of path P(z, w ) can always be normalized to contain only 
non-negative powers of z and w. 





Figure 1: The graphs corresponding to orbifolds a) C 3 /Z 3 x Z 2 , b) the same graph with a different 
choice of path j w , c) C/Z 2 x Z 2 an orbifold of the conifold 

For example, figure [Ija produces the Kasteleyn matrix and determinant 
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(3.4) 
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while the choice of path in figure [l]b (which is chosen so that all crossings occur with positive 
orientation) gives 
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P(z, w) = 1 — 2w + w 2 — 3 zw 2 — 6zw + 3 z 2 w 2 — z 3 w 2 

= w 2 P 0 (3.5) 


The degrees of the monomials of P define points in the lattice Z 2 . The convex hull of 
these points is the Newton polygon of the graph, which provides the link between dimer 
models and toric geometry. 

For a toric three-fold the set of torus actions 


(zi,Z2,Z3) i—► (X ni Zi, X n2 z 2 , X n3 z 3 ) (3-6) 

where A G C*, defines a set of lattice vectors {(ni,n 2 ,ns) G Z 3 } which generate a fan G 
of convex cones. The Calabi-Yau condition is equivalent to the statement that the lattice 
vectors are all coplanar, and their projection onto this hyperplane defines a convex polygon 
in Z 2 , the Newton polygon. 

The generators of the dual cone G define the toric variety algebraically by 


V = Spec(C[A'f nz3 ]) 


(3.7) 


i.e. in terms of the algebraic equations satisfied by the monomials X^X\X^ associated to 
elements (a, b , c) in the dual fan. They correspond to the degeneration locus (“toric skeleton”) 
of the U(l) n torus actions defining the toric variety in the symplectic quotient (linear sigma 
model) construction of V. In string theory language the generators of the dual fan (after 
projection) are the charges (p, q) of 5-branes in the dual brane-web construction of the toric 
variety. See |28|] for further details on (p, q ) 5-brane webs. 

The action of SL( 2, Z) on the torus (i.e. on the periodicity vectors) induces an action 
of SL( 2, Z) on the Newton polygon as a subset of Z 2 . This SL( 2, Z) action is also visible in 
the ( p , q) 5-brane web, where it acts on the charge vector. 


3.1 Examples 

It is a simple exercise to write down some periodic bipartite graphs and compute their 
Kasteleyn determinants as described above. Figure |] shows the important examples of the 
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hexagonal and square graphs (which will be of particular importance to us), and the square- 
octagon graph. Incidentally, these are three of the four different periodic plane tilings with 
even equilateral polygons'. 





Figure 2: The fundamental domain of a) the hexagonal graph corresponding to C 3 , b) the square 
graph corresponding to the conifold C, and c) the square-octagon graph corresponding a toric phase 
of the cone over the Hirzebruch surface Fo . The dashed lines show the periodicity of the graph, 
i.e. how neighbouring domains adjoin. Also marked are an example choice of signs and weights 
associated to the construction of Kasteleyn matrix. 

For the hexagonal graph (which will correspond to C 3 ) the Kasteleyn matrix is 1 x 1, 
and with the choice of signs and paths 7 indicated in figure we find 

P(z, w) = det K(z, w) = 1 — z — w (3.8) 

The corresponding Newton polygon is also shown in figure |3]a together with the absolute 
value of the coefficients. 

Similarly, for the square graph we find 

P(z, w ) = det K(z,w) = 1 — z — w — zw (3.9) 

which agrees with the Newton polygon of the conifold C. 

For the square-octagon graph the Kasteleyn matrix is 


K(z, w) 


P(z, w ) 


/ 1 1 0 ~l/w\ 

1/z -1 1 0 

Owl 1 
v 1 0 -z 1 ) 

det K(z,w ) = — 1/w — w — 1/z — z — 5 


7 We thank R. Kenyon for a correspondence on this point. 


(3.10) 
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a) 1 





Figure 3: The Newton polygon (solid lines)and dual 5-brane web (dashed lines), and the linear 
sigma model field multiplicities for a) C 3 , b) the (resolved) conifold, c) the two toric phases of the 
cone over Fq. 

as shown in figure |3|c. The polygon is equal to the toric diagram of the cone over the 
Hirzebruch surface Fq, and the multiplicity ‘5’ of the internal point agrees with one of the 
two toric phases of this surface, computed in [|T(J (we will explain in the next section how 
to obtain the other phase). This phase is called in the language of || F 2,0 with one double 
impurity. 

3.2 Orbifolds 

Given the fundamental domain of a graph, a natural operation is to enlarge the domain of 
the graph by taking m x n copies of the fundamental domain. It is easily verified that this 
corresponds to subdividing the lattice of the toric diagram, i.e. to orbifolding the space by 
Z m x Z n , where the orbifold action is generated by 


(zi, Z2, Z3) 1—»■ (A^i, Z2, A 1 zs), A" — 1 

(Zi, Z2, Z3) ^ (Zi,(jJZ2,UJ ^3), UJ rn = 1 


(3.11) 

(3.12) 


However, the multiplicities associated to points in the subdivided diagram will also 
change under orbifolding, and it was previously not known how to compute them in terms 
of the un-orbifolded quantities. This problem is solved by using the Kasteleyn matrix of the 
dimer model. The characteristic polynomial of the graph obtained by taking n x m copies 
of a graph is given by PU 


where \ n = l,u; m = 1. 


n— 1 m— 1 

p n , m (z,w )=n n nxz i/n ,u j ™ i/m ) 

i=0 j=0 


(3.13) 
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Notice that in order to use this formula one does not need to know the structure of 
the graph, only its’ characteristic polynomial, i.e. toric data and multiplicities. It can thus 
be applied immediately to compute the multiplicities associated to any orbifold of any non¬ 
compact toric Calabi-Yau for which the multiplicities are known. However in cases where 
the adjacency matrix can be written down for the orbifold in generality (such as for C 3 or 
the conifold, see section [|), it is computationally easier to simply take the determinant of 
this matrix). 

Here are some examples that may be easily calculated using this formula. They agree 
with examples calculated in [19]. |1^, |16[, who also observed the pattern of binomial coefficients 
along the edges of the polygon. This structure is now apparent from formula ( 3.13 ), where 
it arises from the expansion of the product. 


C 3 /(Z 3 x Z 3 ) : 


1-3 3-1 
-3 -21 -3 
3 -3 
-1 


C 3 /(Z 4 x Z 4 ) 


1-4 6-4 1 

-4 -124 -124 -4 
6 -124 6 

-4 -4 
1 


C 3 /(Z 2 x Z 4 ) 


1-2 1 
-4 -12 
6 -2 
-4 
1 


C 3 /(Z 5 x Z 5 ) 


1 

-5 

10 

-10 5 

-5 

-605 

-1905 

-605 -5 

10 

-1905 

1905 

-10 

-10 

-605 

-10 


5 

-5 



-1 





C / (Z 4 x Z 2 ) : 


" 1 -11 
-2 -2 
1 -1 


C/ (Z 3 x Z 3 ) : 


1-3 3-1 

-3 -105 -105 -3 
3 -105 105 -3 
-1 -3 -3 -1 


The above formula may be extended to the case of diagonal Z„ orbifolds by projection of 
the terms appearing in the double product. Faithful actions of the diagonal Z n are generated 
by actions of the form ( 1 , a, b — n — a — 1) (up to permutation), and we find for C 3 


n—1 

P n ( Z , W) = P(A“ W a/n ^ (n - a)/n , A {a+l)i w (a+l)/n z (n-a-l)/n^ ( 3 . 14 ) 

»=0 

where the weights of w and z are related to the coordinates of the corner points that remain 
after orbifolding, see figure [|. 

For example (suppressing minus signs): 
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• The diagonal orbifold C 3 /Z 3 with action (1,1,1), also known as the cone over dP 0 , 
which agrees with 0, 


1 


3 1 
1 


(3.15) 


• The orbifold C 3 /Zg with action (1, 6,1) 

"l 0 0 
0 0 0 
0 0 0 
0 2 0 
0 16 0 
0 20 0 
0 8 1 
0 1 0 


(3.16) 


which is the toric phase of Y' 4,4 |IiJ. Here, and in the following, 0’s indicate points on 
the lattice which are not included in the toric diagram. 


The orbifold C 3 /Z 7 with action (1, 3, 3) 


29 


1 0 0 0 0 
0 7 0 0 0 
0 0 14 0 0 
0 0 0 7 1 
0 0 0 1 0 


(3.17) 


• The other distinct C 3 /Z 7 orbifold with action (1,2,4) |29J 

"l 0 0 0 0 o" 

007000 

000771 

000010 


(3.18) 


The results here are consistent with the analytic formula that was computed for c, the total 
multiplicity of the toric diagrams ( |3.3|) . for an orbifold C 3 /Z n with an action (1,1, —2), | |27| . 
For given n the multiplicity is given by the Fibonacci numbers: 


c = 


'l + y/b 


1 - Vlf 


+ 2. 


(3.19) 
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+ \ 

I \ 

' 

* • » s 



( 0 , 0 ) 

Figure 4: The toric diagram of the diagonal 
orbifold of C 3 /Z n with action (1 ,a,b). 


In terms of the graph, different orbifold ac¬ 
tions correspond to different choices of period¬ 
icity of adjacent copies of the n x m unit cell 
of the orbifold graph, i.e. adjoining neighbour¬ 
ing unit cells with an offset 8 . For example, in 
figure [5] are shown two possible ways of adjoin¬ 
ing neighbouring fundamental domains of the 
Z 2 orbifold of the conifold. These correspond 
to the two different actions of the Z 2 : figure 
| 0 a corresponds to the action ( 1 , 1 , 0 , 0 ) acting 
on the four fields ( w , x, y , z) of the conifold de¬ 
fined by the equation wx = yz, and figure 0]b 
gives the diagonal orbifold action by ( 1 , 1 , 1 , 1 ). 

Computing the Kasteleyn matrix we find 
for figure ||a 

*<*■»> = <»•*» 

P(z,w) = det K(z,w) 

= 1 - w - 2z - 2 wz + z 2 - wz 2 


which agrees with the example computed already in ( |3.14|) , and for figure gb 


8 This looks very similar to the brane box constructions 
similarity; details will appear in the subsequent publication 



Indeed, there is a deeper reason for this 




Figure 5: The two orbifolds of the conifold by Z 2 correspond to different ways adjacent copies of 
the 2 x 1 fundamental domain may be adjoined. 
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(3.21) 


K(z, w ) 
P(z , w) 


/ — I + I/ 2 ; — 1 — l/w\ 

V -1-w 1 - z ) 

det K(z , w) 

— 1/w — w + 1/z + z — 4 


This phase in the language of || is l' 2,0 with two single impurities and can indeed be derived 
as the Z 2 orbifold of the conifold |]3T|. Comparing to (|3.10|) we see that this is another phase 
of the cone over F 0 . The two phases are related by Seiberg duality, as we will explain in 
followup work [^4[ . 

It is often the case that orbifolds of a space appear among the set of phases of larger 
spaces. For example, taking the x Z 3 orbifold of phase ‘5’ of F 0 gives the following 
multiplicities 


-1 

-1 -15 -75 -140 -75 -15 -1 
-1 

and the corresponding orbifold of the phase ’4’ gives 

-1 

-1 -12 -48 -76 -48 -12 -1 
-1 

which both appear among the 18 phases of H 6,0 , listed in appendix |A|. 


(3.22) 


(3.23) 


4. Block determinant formulae for orbifolds 

Since the action of orbifolding a space has a simple realization on the graph, it is often 
possible to write down the general form of the Kasteleyn matrix for arbitrary orbifolds, in 
the form of a block matrix. The blocks are related to the adjacency matrix of the original 
fundamental domain (with a suitable choice of vertex labelling), and the arrangement of 
blocks corresponds to the interconnections between copies of the fundamental domain. Here 
we present the results for the orbifolds of C 3 and of the conifold. Generalizations to diagonal 
orbifold action, and to orbifolds of other spaces are also possible. 

We make use of the block determinant formula for a block matrix M 

m= {cd) ^ 

where A, D are square matrices and A is invertible. Then 
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det(M) = det(.A) det(£> — CA 1 B ) 


(4.2) 



4.1 Orbifolds of C 3 

With edge weights labelled according to 
figure |j. it is easy to verify that the ad¬ 
jacency matrix of the orbifold C 3 /(Z m x 
b„ Z n ) takes the following nm x nm block 
matrix form: 


Figure 6: The edge weights assigned to the (*,j)th 
copy of a fundamental domain. The adjacency ma¬ 
trix for arbitrary multiples of the fundamental do¬ 
main can be written down in block matrix form. 


K(z, w) 


where 


( A 1 Biw 0 . 0 ^ 

0 A 2 B 2 0 0 

0 ... A m _\ B m _i 

\ B m 0 . 0 A m J 

(4.3) 


Ai 


Bi 


( an bn 0 . 0 \ 

0 a i2 b i2 0 ■ ■ • 0 

• ®i(n— 1 ) 1 ) 

\ d . d ^in J 

(cn 0 . 0 \ 

0 Cj 2 0 • ■ ■ 0 


\ o . 0 c in ) 


(4.4) 


Using the block determinant formula (|4.2| ) the evaluation of the determinant can be 
reduced to 


m— 1 


m— 1 


det K = ll det Ai det(kL m - (-1 ) m B m J] A; ] B,) (4.5) 

2—1 2—1 

The following sign choice satisfies the parity rules (IT) and gives the canonical signs for 
the terms appearing in the characteristic polynomial (|3.2|): 


f * bij * 1) C. 




-1 


(4.6) 
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Applying this substitution to K and taking the determinant is a computationally efficient 
way to compute the multiplicities to high rank of the orbifold group. Noting that after the 
substitution (|4.6|) the matrices A, B commute, the determinant (f4.5|) reduces to 


P(z, w) = det(A m - (-1 ) m B m w) (4.7) 

which is easily evaluated. For a visually appealing example, see figure jTj. 

As we discuss in the next section, if the edge weights are kept distinct when taking the 
determinant, one may obtain the phases and multiplicities of toric subdiagrams by Higgsing 
(removing points from the diagram). 


4.2 Orbifolds of the conifold 


where 


K(z, w) = 


(A x z 0 ••• 0 B\ \ 

b 2 a 2 0 0 

0 B 3 A 3 0 0 

Vo 0 B m A m J 


A, = 


B, = 


(Cji b n w 0. 0 \ 

0 Ci 2 bi 2 0 • ■ ■ 0 

■ Cj(n—1) ^i(n—1) 

\ bi n 0 . 0 Ci n J 

(dn aniu 0 . 0 \ 

0 di 2 cii2 0 • • • 0 

• * * • di(n— 1 ) 1 ) 

\ ^in 0 . 0 di n J 


An appropriate sign choice is 


d-ij > 1 1 bij * 1 ? C{j * 1 ? dij 1 


(4.8) 


(4.9) 


(4.10) 


As in the previous section, the determinant of the block matrix can be simplified with this 
sign choice, giving 


P(z,w) = det (A m - (-1 ) m B m z) 


(4.11) 
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5. Removing points from the toric diagram (partial resolution) 


Given a toric diagram, one may obtain other toric diagrams as subdiagrams by removing 
points from the corners (this preserves convexity of the polygon). Given the interpretation 
of P as counting of paths on the graph, in order to remove a monomial from P we must 
cut all paths in the graph with the specified height changes. This can be done by removing 
certain edges from the graph, i.e. setting their weights to zero. There is often a choice of 
which edge(s) to cut to implement removal of a given set of points from the toric diagram, 
and one finds that the set of multiplicities of the points in the resulting diagram enumerate 


the possible phases of the toric space that may be obtained by Higgsing, as studied in [16 


Certain graphs, such as the hexagonal and square graphs, are universal in the sense 
that any other toric diagram may be obtained from them by taking a suitably large orbifold 
by Z m x Z n and then removing points from (Higgsing) the resulting toric diagram. This 
provides a systematic computational procedure for investigating the phase structure of the 
associated gauge theories. In practise, in order to simplify the calculation it is desirable to 
choose an minimal embedding; i.e. an embedding into a larger graph with a minimal number 
of additional points. The following algorithm may be used: 

The adjacency matrix of the embedding graph is written with all edge weights distinct 
and its’ determinant is computed. Then evaluation of the phases of the toric diagram may 
be achieved by solving the set of simultaneous equations specifying that the coefficients in 
the determinant of the points to be removed from the toric diagram all vanish, while the 
coefficients of the remaining points do not. Once the set of solutions has been found, the 
remaining edge weights may be set to a suitable choice of ±1 to recover the set of phases 
and multiplicities of the target diagram. 

A refined version of this algorithm is as follows: 


1 . Compute the characteristic polynomial of the ambient toric diagram with all edge 
weights parameterized distinctly (as in Figure |6|). 

2. Identify the set of vertices to be removed from the ambient toric diagram, ordered so 
that the polygon remains convex after each removal. These are by assumption corner 
points of the toric diagram at each stage of the removal. 

3. Take the complement of this set, i.e. identify the points in the diagram that we do not 
wish to remove. For each point (corresponding to a monomial in the determinant), 
identify the edge weights in the coefficient of this monomial that are common to all 
terms, i.e. which factor out. We cannot set any of these variables to zero, or an 
unwanted extra point will be removed from the toric diagram. 
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4. For each vertex in the list to be removed, identify the edge weights in the coefficient 
of the determinant that are common to all terms, i.e. which factor out. Discard any 
edge weights that are on the list computed in step 


5. For each variable in this list, set it to zero to remove the point from the toric diagram, 
and recurse to step f| with the new set of corners of the sub-diagram. 


6 . If all possible cuts have been made (i.e. no more variables remain to be zeroed), and 
we have removed all of the points in the target list, we have found a solution. The 
edges that were cut to obtain this graph are known, and the multiplicities of this phase 
may be read off by evaluating the determinant against the variables set to zero and 
the canonical sign assignment of the initial graph (e.g. (|4.6|) for C 3 ). 


The limitations of this algorithm are discussed in the following section. One advantage of 
the algorithm compared to methods discussed in previous sections is that it operates on the 
graph, not the determinant (which contains much less information). Thus, solutions come 
not just with the toric multiplicities but the subgraph of the initial graph which produces 
them. 

Note that sometimes it is possible to obtain the desired shape of the toric diagram such 
that one or more of the corner points have multiplicity greater than one. This typically 
corresponds to an inconsistent quiver gauge theory. However, in all examples we have en¬ 
countered, it is always possible to perform additional edge cuts to reduce this multiplicity 
to 1 (this is taken care of by the form of the algorithm presented above). 

Using this method we are able to recover all previously known results on multiplicities 
of toric diagrams (low-rank orbifolds of C 3 and the conifold; cones over F 0 , dPo,i, 2 , 3 , the 
toric spaces Y p ’ q and X p,q , and some other miscellaneous examples), as well as generating 
many new examples which have not previously been studied. Some simple new examples are 
presented here: 


• The orbifold C 3 /(Z 2 x Z 4 ) with point (0,4) removed. 


1-2 1 
-3 -9 
3 -1 
-1 


(5.1) 


• The orbifold C 3 /(Z 4 x Z 4 ) with points (0,4), (4,0) removed. There are 6 phases: 
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1-3 3-1 

-3 -65 -43 -1 
3 -51 2 
-1 -1 

" 1 -3 3-1 

-3 -65 -51 -1 
3 -43 2 
-1 -1 

It is interesting to observe the symmetry in the above set of phases; two pairs of phases are 
interchanged by reflection, and the other two phases are symmetric. This structure should 
also be observed in the phases of the dual quiver gauge theory (the phases of the toric space 
F 6 ' 0 listed in the appendix exhibit a similar structure). 

Other toric diagrams may be computed similarly (Mathematica code available upon 
request). 

6. Discussion 

We have proposed the existence of a duality between quiver gauge theories and dimer mod¬ 
els, and developed dimer model techniques for enumerating the toric phases of the gauge 
theories, and for computing their linear sigma model field multiplicities. We obtained general 
formulae for the multiplicities of arbitrary orbifold theories, and discussed how to obtain the 
multiplicities of non-orbifold phases by Higgsing a larger toric diagram. These techniques 
reproduce all known results from existing quiver gauge theories, and give many predictions 
for quiver theories not previously examined. 

As currently understood, the Higgsing algorithm contains two computational limitations: 

1. The need to take the determinant of the orbifold adjacency matrix with all edge weights 
distinct. Even making use of the block structure of the orbifold matrix as in section 
[i], it becomes prohibitively difficult to calculate the determinant for higher orbifold 
rank (e.g. beyond C 3 /(Z 6 x Z 6 ) or so, i.e. a 36 x 36 matrix with 3 x 36 = 108 
distinct edge weights). With a better understanding of the relationship between a 
given toric diagram and the properties of the dimer model(s) that generate it, it may 
be possible to identify the edges to be removed without needing to solve the constraint 
equations, in which case it would not be necessary to use distinct edge weights during 
the computation. As figure [7| suggests, this would allow the procedure to be extended 
to ambient spaces with much higher orbifold ranks.) 


1-3 3-1 


1-3 3-1 

-3 -71 -53 -1 


-3 -81 -59 -1 

3 -45 2 


3 -59 2 

-1 -1 


-1 -1 

" 1 -3 3 -l" 


" 1 -3 3 -l" 

-3 -71 -45 -1 


-3 -73 -43 -1 

3 -53 2 


3 -43 2 

-1 -1 


-1 -1 
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2. The brute-force enumeration of solutions, which does not account for the large amount 
of symmetry among the set of solutions. For example, a set of removed edges may 
often be translated around the graph without changing the characteristic polynomial. 
In practice this means that the running time of the algorithm grows quickly with the 
number of points to be removed from the diagram, because of a large overcounting 
of solutions. It would be useful to precisely characterize and remove the redundancy 
among solutions. 


Typically one finds that the algorithm produces many of the phases very quickly, with 
the exhaustive search for all phases taking longer. A shortcut would be to take these 
first solutions and use Seiberg duality to fill out the set of phases. The interpretation 
of Seiberg duality in dimer language will be discussed in |24| . 


The reason for the correspondence between quiver gauge theories and dimer models, and 
the underlying string theoretical interpretation will be discussed in forthcoming work [ p4| . 
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A. Toric phases of Y 6,0 


The following are the multiplicities of 17 of the 18 phases [^3j] of the toric space F 6,0 (sup¬ 
pressing the minus signs from ( |3.2| )), and the edges in the graph of the Z 2 x Z 6 orbifold of the 
conifold (using the labeling of figure |B| and the block matrix described in section [4.2p that 
may be cut to obtain this phase. Many other equivalent choices of edges are possible. The 
“missing” 18th phase was identified in [|3T| and shown to have identical toric multiplicities 
to one of the remaining 17. 


1 

1 14 61 102 65 15 1 
1 

1 

1 15 65 102 61 14 1 
1 

1 

1 14 63 106 67 15 1 
1 

1 

1 15 67 106 63 14 1 
1 

1 

1 12 48 76 48 12 1 
1 

1 

1 13 54 86 54 13 1 
1 

1 

1 13 54 87 54 13 1 
1 


<711, <721, <731, a41, a51, a61, cl2, 522, c42, 532, c62, 552 


<711, <721, <731, a41, a51, a61, c22, 512, c32, 542, c62, 552 


<711, <721, <741, a31, a51, a61, cl2, 522, c42, 532, c62, 552 


<711, <721, <741, a31, a51, a61, c22, 512, c32, 542, c62, 552 


<711, <721, <731, a41, a51, a61, c32, 512, c42, 522, c52, 562 


<711, <721, <731, a41, a51, a61, cl2, 522, c32, 552, c42, 562 


<711, <721, <731, a41, a51, a61, cl2, 522, c32, 542, c52, 562 
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1 

1 13 55 89 55 13 1 
1 

1 

1 14 62 104 62 14 1 


1 14 63 102 63 14 1 
1 

1 

1 14 63 103 63 14 1 
1 

1 

1 14 64 109 64 14 1 
1 

1 

1 15 63 100 63 15 1 
1 

1 

1 15 65 105 65 15 1 
1 

1 

1 15 75 140 75 15 1 
1 

1 

1 16 78 132 78 16 1 
1 

1 

1 18 79 126 79 18 1 
1 


dll, d21, <731, a41, a51, a61, c22, 612, c42, 632, c52, 662 


dll, d21, d31, a41, a51, a61, c22, 612, c42, 632, c62, 652 

dll, d21, d41, a31, a51, a61, c32, 612, c42, 622, c62, 652 

dll, d21, d41, a31, a51, a61, cl2, 622, c32, 652, c42, 662 

dll, d21, d41, a31, a51, a61, cl2, 622, c32, 642, c52, 662 


dll, d21, d31, a41, a51, a61, cl2, 622, c52, 632, c62, 642 

dll, d21, d31, a41, a51, a61, cl2, 622, c32, 642, c62, 652 

dll, d31, d51, a21, a41, a61, c22, 612, c42, 632, c62, 652 

dll, d21, d41, a31, a51, a61, cl2, 622, c32, 652, c42, 662 


dll, d21, d31, a41, a51, a61, cl2, 622, c32, 642, c62, 652 

(A.l) 


22 



References 


[1] 

[ 2 ] 

[3] 

[4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 

[ 11 ] 

[ 12 ] 

[13] 

[14] 

[15] 

[16] 


A. Okounkov, N. Reshetikhin, and C. Vafa, Quantum Calabi-Yau and classical crystals, 

|hep-th/030920£ . 

A. Iqbal, N. Nekrasov, A. Okounkov, and C. Vafa, Quantum foam and topological strings, 

|hep-th/0312022 . 

J. Stienstra, Mahler measure variations, Eisenstein series and instanton expansions, 

frath. NT/0502193 . 


J. Stienstra, Mahler measure, Eisenstein series and dimers, nath.NT/0502197 . 

S. Franco, A. Hanany, and Y.-H. He, A trio of dualities: Walls, trees and cascades , Fortsch. 
Phys. 52 (2004) 540-547, ||hep-th/0312222|] . 

S. Benvenuti, A. Hanany, and P. Kazakopoulos, The toric phases of the Y(p, q ) quivers, 

|hep-th/0412279| . 

D. R. Morrison and M. R. Plesser, Non-spherical horizons. I, Adv. Theor. Math. Phys. 3 
(1999) 1-81, jhep-th/9810201 1. 

C. Beasley, B. R. Greene, C. I. Lazaroiu, and M. R. Plesser, D3-branes on partial resolutions 
of Abelian quotient singularities of Calabi-Yau threefolds, Nucl. Phys. B566 (2000) 599-640, 

Hhep-th/9907186|1 . 

B. Feng, A. Hanany, and Y.-H. He, D-brane gauge theories from toric singularities and toric 
duality, Nucl. Phys. B595 (2001) 165-200, | hep-th/0003085 |. 

C. E. Beasley and M. R. Plesser, Toric duality is Seiberg duality, JEEP 12 (2001) 001, 
Hhep-th/0109053| 1. 

B. Feng, A. Hanany, Y.-H. He, and A. M. Uranga, Toric duality as Seiberg duality and brane 
diamonds, JHEP 12 (2001) 035, |hep-th/0109063 ]. 

F. Cachazo, B. Fiol, K. A. Intriligator, S. Katz, and C. Vafa, A geometric unification of 
dualities, Nucl. Phys. B628 (2002) 3-78, ||iep-th7oi l002S|. 

A. Hanany and A. Iqbal, Quiver theories from D6-branes via mirror symmetry, JHEP 04 
(2002) 009, [|hep-th/0108137|1 . 

C. P. Herzog and J. Walcher, Dibaryons from exceptional collections, JHEP 09 (2003) 060, 

Hhep-th/0306298|1 . 

P. S. Aspinwall and S. Katz, Computation of superpotentials for D-branes, [hep-th/0412209| . 

B. Feng, S. Franco, A. Hanany, and Y.-H. He, Symmetries of toric duality, JHEP 12 (2002) 
076, Hhep-th/0205144|1 . 


23 





































[17] T. Muto and T. Tani, Stability of quiver representations and topology change, JHEP 09 
(2001) 008, [|hep-th/0 107217|] . 

[18] B. Feng, A. Hanany, and Y.-H. He, Phase structure of D-brane gauge theories and toric 
duality, JHEP 08 (2001) 040, [[hep-th/ 0 104259 1. 

[19] A. Dunn and A. Hanany. Work in progress. 

[20] R. Kenyon, A. Okounkov, and S. Sheffield, Dimers and amoebae, tnath-ph/ 0311005[ 

[21] T. Sarkar, D-brane gauge theories from toric singularities of the form C 3 /r and C 4 /T, Nucl. 

Phys. B595 (2001) 201-224, [|hep-th/0005166|1 . 

[22] S. Benvenuti, S. Franco, A. Hanany, D. Martelli, and J. Sparks, An infinite family of 
superconformal quiver gauge theories with Sasaki-Einstein duals, li ep-th/041 1 264 . 

[23] A. Hanany, P. Kazakopoulos and B. Wecht, A new infinite class of quiver gauge theories 

|hep-th/05031 77. 

[24] S. Franco, A. Hanany, K. D. Kennaway, D. Vegh and B. Wecht, Brane dimers and quiver 
gauge theories , hep-th/0504110. 


[25] R. Kenyon, An introduction to the dimer model, nath. C0/0310326 . 

[26] P. Kasteleyn, Graph theory and crystal physics, in Graph theory and theoretical physics, 
pp. 43-110. Academic Press, London, 1967. 


[27] T. Muto, D-geometric structure of orbifolds, [hep-th/0206012 


[28] O. Aharony, A. Hanany, and B. Kol, Webs of ( p, q) 5-branes, five dimensional field theories 
and grid diagrams, JHEP 01 (1998) 002, [hep-th/9710116]. 

[29] Both of the Z 7 toric diagrams were checked by Yang-Hui He in an unpublished work. 

[30] A. Hanany and A. M. Uranga, Brane boxes and branes on singularities, JHEP 05 (1998) 013, 
[|hep-th/9805139|l . 

[31] S. Franco, A. Hanany, and A. M. Uranga, Multi-flux warped throats and cascading gauge 


theories, [hep-th/05021 13 


24 



























s 




S 

a ^ 


... 


Figure 7 : Multiplicities for the orbifold C 3 /(Zio x Z50). 
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